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t— I ■ Abstract 

If} \ We study the production of the D and J/ tp mesons from deconfined quark mat- 

CTN ■ ter at CERN SPS energy. Using the MICOR microscopical coalescence model we 
determine the transverse momentum spectra of these charm mesons. We predict the 

^) \ slopes of their transverse momentum spectra in Pb+Pb collision at 158 GeV/nucleon 

' beam energy. 

o 1 

^ ' During the last decade charm hadron production was investigated with great inter- 

est at the CERN Super Proton Synchrotron (SPS), because the theoretically predicted 
^h. anomalous suppression of the J/ tp particle could indicate the appearance of colour decon- 

finement and the formation of the quark-gluon plasma (QGP) phase M. An anomalous 



suppression has been detected in Pb+Pb collisions at 158 GeV/nucleon bombarding en- 
ergy |H| and this result is considered as the strongest evidence of QGP formation (for an 



overview see Ref. M). However, the unambiguous confirmation of the QGP formation re- 



quires more charm related data, e.g. D-meson spectra and absolute yields. This provides 
motivation to measure open charm production in Pb+Pb collision at the CERN SPS ||. 
In parallel, new theoretical efforts work to determine charm hadron yields and ratios in 
this energy range using statistical |fl , thermal |6j and non-equilibrium J7| models of charm 
hadron production. 

Since momentum spectra could reveal some interesting details about the formation 
mechanisms of charm hadrons and thus the properties of the early (supposedly deconfined) 
phase, our aim is to determine the transverse momentum spectra of the D and J/ip mesons 
in heavy ion collisions. We assume that hadrons are produced directly from a deconfined 
state by quark coalescence. We neglect hadronic rescatterings. This hadronization sce- 
nario is modelled using Microscopical COalescence Rehadronisation (MICOR) || 10 



which has been successfully applied to reproduce the measured transverse slopes of the 
0, Q and p mesons p[ . (The spectra of other light and strange hadrons were investigated 
previously in Refs. 00.) The MICOR model is the successor of the ALCOR [11], p 



and the Transchemistry |13[ models. While the ALCOR and Transchemistry models were 
constructed to determine the total number of hadrons produced from quark matter, the 



MICOR model can also determine the full momentum spectra of the produced hadrons. 
We therefore include charm quarks into the MICOR model and calculate the spectra of 
directly produced charm mesons. 

In the MICOR model we start with a deconfined initial phase as motivated by lattice 
QCD calculations. (Note that phase differs from the idealized, asymptotically free quark- 
gluon plasma [00 The lattice results for the QCD equation of state (with Nf = 0,2,4 
flavour) were parametrized phenomenologically as a massive quark matter JT5fl . This 
parametrization reproduces the lattice QCD results regarding the pressure and energy 
density in the temperature region 1 < T/T c < 3 |T5] and thus it can be used as a suitable 
description of hadronization. The massive gluons (whose number is suppressed relative to 
the quarks and antiquarks) create an attractive effective force between the colored quarks 
and anti-quarks which drives the hadronization. According to our assumption, the quark 
matter here is a fully thermalized state, and in the hadronization stage is characterized by 
a hadronization temperature Thad, a transverse flow and a Bjorken-scaled longitudinal 
flow. 

The coalescence of the massive quarks produces a hadron resonance gas, which is out 
of equilibrium || . This out of equilibrium character will not change during the resonance 
decay. While the final state hadronic spectra can be fitted by exponential functions on 
transverse momenta, this does not mean that the obtained slope parameters characterize 
equilibrated hadronic gas states. 

During hadronization meson-like objects are formed via quark- ant iquark coalescence 
(see Ref. for the coalescence cross section), while baryon-like objects are formed in two 
steps: the formation of diquarks from quarks are followed by the coalescence of diquarks 
and quarks into baryons. We assume that color-neutral prehadrons escape the deconfined 
state without desintegration and that they become part of the produced excited hadronic 
gas. These prehadrons are off-shell, however their off-shell masses are very close to the 
well-known bare hadron masses. This minor difference is corrected in a last step of the 
hadronization, when a prehadron absorbs or radiates the appropriate (small) energy to 
become on-shell. We assume that it conserves its velocity during this process. More 
details about the physics of the MICOR model can be found in Ref. ||. 

The initial condition of the MICOR hadronization model consists of the hadroniza- 
tion temperature, Thad-, and the initial transverse flow, v^, of the deconfined phase. The 
analysis of the measured transverse momentum spectra of 0, Q and p particles was per- 
formed in Ref. || and fits to these experiment data favor initial values of Thad — 175 ±15 
MeV and of v% = 0.46 ± 0.05 in the Pb+Pb collision. We will use the mean values of 
these parameters to calculate charm meson production. The theoretical uncertainty of 
our calculation will be estimated through the uncertainty of these initial parameters. We 
consider a large enough longitudinal extension for the deconfined matter: r] = ±2.2, 
where r] is the space-time rapidity (this parameter will not influence the transverse mo- 
mentum spectra at y = 0). The constituent quark masses are chosen to be rriQ = 310 
MeV, m s = 430 MeV, m c = 1540 MeV. 

We calculate the transverse momentum spectra of D and J /if) mesons in Pb±Pb 
collision at Thad — 175 MeV and = 0.46. The D meson is produced from the coalescence 
of a charm quark and a Q antiquark, while the J /if) is produced as a CC state. Figure 
1 displays the obtained spectra with arbitrary normalization. The full line shows the 



MICOR result and the dashed line indicates the best fit with the parametrization 



dN / mTdrriT = C ■ m^ 2 exp(— /T e f f) (1) 

in the momentum regions 0.3 GeV < — < 2.2 GeV . In this region we obtain 
Tf f b f +Pb (D) = 259 MeV and Tf f b f +Pb (J/^) = 315 MeV for the effective slope parameters. 
Figure 1 shows how the non-thermal momentum distribution produced in the coalescence 
process (full curve) mimics the thermalized distribution of eq. ([I]) (dashed curve) in large 
momentum region. 

The NA38 Collaboration has measured the transverse spectra of the J/ip in S+U 
collision at 200 GeV/nucleon bombarding energy [16|]. To apply our model, we assume 
that the massive quark-antiquark matter was already produced in S+U. In Figure 1 the 



experimental data from Ref. [16] are marked by full dots (at arbitrary normalization). The 
full curve is the MICOR calculation using T had =175 MeV and a reduced transverse flow, 

= 0.33. This reduced value for the transverse flow can be understood, since the smaller 
project radius of S+U collision yields a smaller transverse size and shorter expansion 
time for the deconfined phase as compared to that in the Pb+Pb collision. Using the 
parametrization of eq. ([I]) in the momentum region 0.3 GeV < nix — mj/^ < 2.2 GeV we 
obtained T^ u (J/ijj) = 234 MeV (dashed line) thus recovered the experimental value, 
T^ u (J/ijj) = 234 + 3 MeV, obtained in the highest transverse energy bin (see Ref. JIB]). 
(The published slope value was obtained with the fit dN / mTdrriT = C ■ rrirp ■ ^(my/T), 
which is approximately the same one as eq. ([!]) in the temperature region T « 200 — 300 
MeV for particles heavier than the proton, see Ref. ||.) For the D-meson we obtain 
Teff U (D) = 211 MeV. Figure 1 shows the MICOR result (full line), while the fit from 
eq.(|l|) is indicated by dashed line. 

We now investigate the dependence of the final transverse slope on the hadronization 
temperature and transverse velocity. First we fix the hadronization temperature at T^ad = 
175 MeV and vary the transverse flow parameter, v^. Figure 2 shows our predictions on 
the effective slopes for the J / ip (upper curve) and for the D meson (lower curve) spectra. 
The open squares show the MICOR results in the Pb+Pb collision at 158 AGeV energy 
at Vj> = 0.46. The uncertainty in the transverse flow, Avj, = ±0.05, yields an errorbar on 
the theoretical predictions: T Pb f +Pb (D) = 259 ± 19 MeV, T Pb f +Pb (J/ip) = 315 ± 35 MeV. 
The NA38 result on J/ip in S+U collision is indicated by the star and open circles show 
the MICOR predictions at Vj, = 0.33 . Assuming a 10 % uncertainty in the transverse 
flow, namely At>y = ±0.03, we obtain the theoretical predictions: T^ U (D) = 211 ± 9 
MeV, T^ff 11 (J / 4>) = 234 ± 13 MeV. These uncertainties show the strong dependence of 
the charm hadron slopes on the transverse flow of the intermediate deconfined phase. As 
new experimental data will appear, analysis of Figure 2 can be used to determine the 
appropriate parameter for the MICOR model at the Thad = 175 MeV hadronization 
temperature. 

We repeat this calculation also varying the hadronization temperature, T^ad-, with a 
AThad = 15 MeV. The obtained theoretical uncertainty remains in the same scale, ±15 
MeV. 

In Figure 3, we display the MICOR results for the T e ff slopes of the D and J / ip mesons 
produced in Pb+Pb collisions (they are marked with open squares). We also display the 
measured and the calculated slope parameters for the p, cf> and Q M. Experimental data 



on the slope parameter of other particles are from Refs. [17, 18, [T^, ^(J. (All data was 
parametrized by eq. (|l|) in the momentum region — rrii > 0.3 GeV. ) 

A similar plot can be created for the S+A collisions at 200 GeV/nucleon project energy. 



Figure 4. includes results for S+U [fLq[ , S+W |22] and S+Au |23| collisions. The open 
squares show the MICOR results for particles p, (ft, Q, D and J/ip. Here the J/ip data 
determines the transverse flow parameter, v ^ = 0.33, which reproduce the slope for the 4> 
meson. 

Figures 3 and 4 support the assumption that the secondary hadronic interactions only 
slightly modify the transverse momentum spectra at higher momenta, my — > 0.3. 
They reveal that the separation of the p + ,p~, A, A, S~ and H + from the weakly interacting 
and Q is small. Even more, the slopes of the p + and p~ are very close to that of the 
weakly interacting meson. Thus we do not expect considerable modification of the much 
heavier charm mesons transverse spectra. 

Our result sides with the analysis of hydrodynamical expansion versus hadronic cas- 
cade evolution performed by URQMD [p5| , where enhanced initial transverse flow was 



paired with a reduced influence of hadronic cascading in the transverse momentum spec- 
tra of heavier particle. In this manner, our result does not contradict to recent finding of 
Lin et al. ]2E|], where relatively large modification was found from "hadron+D" collisions 
of the D-meson spectra in a longitudinally expanding system. Allowing a relatively large 
transverse flow, e.g. = 0.46 as in MICOR, the lifetime of the fireball will be shortened 
and the efficiency of the hadronic rescattering will be reduced. 

In this paper, we predict the transverse slope of the D and J/ip mesons in Pb+Pb 
collision at CERN SPS energy. Using the microscopical MICOR model we could follow the 
hadronization of the quark matter that we assume is produced in this heavy ion collision, 
and calculate the transverse momentum distributions. We obtained T Pb j~ Pb (D) = 259±19 
MeV and T Pb j~ Pb (J/ip) = 315 ± 35 MeV. We compare our model with experimental data 
on J I ip production in S+U collision and demonstrate that we can reproduce the measured 
spectra with a reasonably reduced transverse flow of the quark matter at a hadronization 
temperature of Thad = 175 MeV. We predicted the slope parameter of the D-meson to 
be T^jf U (D) = 211 ± 9 MeV. We urge future experiments to measure this quantity to 
provide further tests of this model. 
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Fig. 1. The transverse momentum spectra of the D and J/ip meson calculated by 
the MICOR model with parameters T had = 175 MeV and «§. = 0.46 for Pb+Pb 
collision at 158 AGeV and with v% = 0.33 for S+U collision at 200 AGeV (solid 
lines). The dashed lines are the fits of eq.(|l|) in the transverse momentum region 
0.3 < itlt — rrii < 2.2 GeV. The full dots are the data from the NA38 Collaboration 
for S+U collision flTC 
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Fig. 2. Prediction from the MICOR model for the T e ff slopes of the D and J/i/j 
transverse momentum spectra as the function of initial flow at = 175 MeV. 
The open squares indicate the MICOR prediction with = 0.46 for Pb+Pb colli- 
sions at 158 AGeV and the open circles indicate the prediction with v j< = 0.33 for 
the S+U collisions at 200 AGeV. The errorbars shows the influence of 10 % uncer- 
tainty in the transverse flow. The star represents the NA38 experimental data for 



the J/ip slope in S+U collisions, see Figure 1 and Ref. [|K 
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□ MICOR model: T H =175 MeV, v T =0.46, (p, 4>, Q, D, J/ty) 

Exp. data in Pb+Pb (158 AGeV): 

■ WA97 - SQM98, JPG25(1 999)1 81 

• NA50 - QM99, NPA661 (1 999)534 

A NA44 data, PRL78(1 997)2080 

fit! dN/m T dm T = C m T 1/2 exp(-m T /T eff ) 
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Fig. 3. Experimental hadronic slopes of the transverse momentum spectra in the 
PbPb collision at 158 AGeV energy from WA97 |T7| (squares), NA50 f0§ (dots) 
and NA44 Collaboration |19j (triangulars) . The results of WA97 and NA50 were 



originally fit by eq. ([[]). We fit the NA44 data on tt ± , K ± , p + and p |T9l in the 



same way in the momentum region tut — rrii > 0.3 GeV. Open squares indicate the 
MICOR results on charm mesons and other hadrons. 
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□ MICOR model: T H =175 MeV, v T =0.33, (p, 4>, Q, D, J/f) 

Exp. data in S+W (200 AGeV): 

■ WA85 - PLB355,B359,B376,B412 

• NA38 - PLB368(1 996)239, PLB423(1 998)207 

A NA35 data in S+Au (200 AGeV) 

fit! dN/m T dm T = C m T 1/2 exp(-m T /T eff ) 
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Fig. 4. Experimental hadronic slopes of the transverse momentum spectra in the 



S+W and S+Au collisions at 200 AGeV energy from WA85 0, g2j (squares), NA38 
H| m (dots) and NA35 Collaboration |3| (triangulars) . The results of WA85 and 



NA38 data were originally fit by eq. ([]]). We fit the NA35 data on 7r~, K° s , p + , A 



and A p3| in the same way in the momentum region — m ; > 0.3 GeV. Open 
squares indicate the MICOR results on charm mesons and other hadrons. 



